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We analyse simulations of turbulent, magnetised molecular cloud cores focussing on the for¬ 
mation of Class 0 stage protostellar discs and the physical conditions in their surroundings. We 
show that for a wide range of initial conditions Keplerian discs are formed in the Class 0 stage 
already. Furthermore, we show that the accretion of mass and angular momentum in the sur¬ 
roundings of protostellar discs occurs in a highly anisotropic manner, by means of a few narrow 
accretion channels. The magnetic field structure in the vicinity of the discs is highly disordered, 
revealing field reversals up to distances of 1000 AU. These findings demonstrate that as soon as 
even mild turbulent motions are included, the classical disc formation scenario of a coherently 
rotating environment and a well-ordered magnetic field breaks down. 


1 Introduction 

In our research we investigate the formation of protostars and their associated protostellar 
disks, the early precursors of planetary systems. During the last decade simulations of 
collapsing molecular cloud cores have revealed the so-called catastrophic magnetic braking 
problem: Magnetic fields are able to transport angular momentum by means of toroidal 
Alfven waves. Modeling the collapse of rotating molecular cloud cores, simulations have 
shown that in the presence of magnetic fields with strengths comparable to observational 
results, the formation of rotationally supported (Keplerian) protostellar disks is largely 
suppressecl31. This is due to the fact that angular momentum is removed very efficiently 
from the interior of the core by the magnetic field. In previous works we could confirm 
this effect for the collapse of massive (100 solar masses), molecular cloud coreP. This 
key result of the suppression of Keplerian disk formation during the earliest stages of star 
formation is in contrast to recent observational results which state that protostellar disks 
should be present already in the Class 0 stagd33 

2 Numerical methods 

The simulations presented here are performed with the hydrodynamics code FLASH, 
which is written in Fortran 9CP! The code solves the 3-dimensional, discretized mag- 
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netohydrodynamical equations on a Cartesian grid. Making use of the adaptive-mesh- 
refinement (AMR) technique, only those regions which are of particular interest for us are 
resolved with the highest possible spatial resolution whereas other regions of minor inter¬ 
est are resolved more coarsely. This significantly reduces the number of calculations to 
be performed and hence the computational time required, thus allowing us to perform the 
simulations over long physical timescales. We also make use of the sink particle routine to 
model the formation of protostarJ^. 

3 Initial conditions 

Observations of the birth places of stars show a wide range of physical quantities, in partic¬ 
ular in their initial mass. As we do not simulate a particular region observed by astronomers 
but rather aim to understand the systematic influence of the initial conditions, we have to 
perform a number of simulations in our work covering a wide range of masses and turbu¬ 
lence strengths. This allows us to draw conclusions about the effect of the initial conditions 
on the formation mechanism of stars. We modeled the collapse of molecular cloud cores 
with masses ranging from about 2 solar masses up to 1000 solar masses. The cores are 
threaded by a strong magnetic held along the z-axis and have an additional supersonic, 
turbulent velocity held as indicated by observations. 

4 Results 

The results of previous simulations described in section 1 show up in case that highly 
idealized initial conditions are used for the simulations. In particular the lack of turbulent 
motions - frequently observed in molecular cloud cores - could have a significant effect on 
the formation of protostellar disks and outflows. For this reason, in our research we here 
focus on the influence of turbulence on the formation of protostellar disks and outflows. 
This work has been performed on JUROPA and other supercomputing facilities. Each 
of the simulations required a computational time of a few 100 000 CPU-hours with a 
simultaneous use of up to 1000 CPUs per simulation. A few hundreds of files were produce 
for each simulation requiring a disk space of a few TB in total. 

4.1 Turbulence-induced disk formation 

We examine our simu lations focusing on the question of how turbulence affects the forma¬ 
tion of Keplerian diskJ33i3 An example result is shown in Fig. 1 showing the protostellar 
disk in a representative run with a molecular cloud core of 100 solar masses. 

As can be seen, in the presence of turbulence rotationally supported disks are re¬ 
obtained again, which is in contrast to the previous simulations with comparable magnetic 
held strengths but no turbulencJ^. This suggests that the efficiency of the magnetic braking, 
which is responsible for removing angular momentum from the midplane, is reduced sig¬ 
nificantly in the turbulent case. Analyzing the surroundings of the disks we can show that 
this indeed the case. The turbulent surroundings of the disk exhibit no coherent rotation 
structure (compare Fig. 1). Therefore, an efficient build-up of a strong toroidal magnetic 
held responsible for angular momentum extraction is hampered. Moreover, the turbulent 
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Figure 1. Protostellar disk seen from top-on. Black dots represent protostars, green arrows the velocity field. 



motions lead to a strongly disordered magnetic field which further reduces the magnetic 
braking efficiency. Since simultaneously the angular momentum inwards transport remains 
high due to the presence of local shear flows in the vicinity of the disks, there is a net in¬ 
wards angular momentum transport towards the center of the disk. The combination of 
these effects results in the observed build-up of Keplerian disks as expected from obser- 
vationJ33 Varying the core masses (2.6 - 1000 solar masses) and the turbulence strengths 
does not change our findings. This clearly demonstrate that the turbulence-induced disk 
formation mechanism works for a wide range of initial conditions. In particular, we could 
show that the formation of Keplerian disks does not require an uniform rotation of the core 
- turbulent motions alone lead to the build-up of Keplerian disks. Moreover, we showed 
that even for subsonic turbulence, which is usually present in low-mass protostellar cores, 
the turbulence-induced formation mechanism still holds. 

In Fig. 2 we show the 3-dimensional structure of the magnetic field and the gas motions 
around a protostellar disk in one of our simulations. The magnetic field lines clearly reveal 
a highly complex structure being far off from well-ordered. Several field reversals up to 
distances of about 1000 AU from the disk center can be found. This indi cates th at the 
recently observed inclination of magnetic fields with respect to the disc axi£3SEl could 
simply be the consequence of a spatially unresolved, highly disordered field structure. 
The anisotropy of the accretion towards the disk is represented by the colored regions 
in the figure. Accretion is dominated by a few narrow accretion channels which carry a 
significant part of the inflowing mass and angular momentum but cover only about 10% of 
the surface area. We emphasize that we find this anisotropic accretion mode in all of our 
simulations. 

5 Conclusions and outlook 

To summarize, our simulations present a richer picture of the process of disk formation, 
wherein turbulence, filamentary accretion streams, and magnetic field reversals guarantee 
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Figure 2. 3-dimensional stnacture of the magnetic field and gas motions around a Keplerian disk in one of our 
simulations. 


that the otherwise overwhelming strength of magnetic braking by ordered fields is signifi¬ 
cantly degraded, allowing Keplerian discs to form. Non-turbulent collapse simulations, on 
the other hand, might significantly overestimate the efficiency of magnetic braking and thus 
underestimate the fraction of Class 0 stage Keplerian discs. We suggest that the anisotropic 
accretion and disordered magnetic held structure found in the environment of protostellar 
discs might set the stage for other mechanisms c ontrib uting to the formation of discs like 
(pseudo-) disc warping or non-ideal MHD effect£5IIIl 

For our future research we intend to study the self-consistent ejection of protostellar 
outhows from the Keplerian disks formed in our turbulence simulations. In order to reach 
this goal we will have to redo some of our simulations with increased spatial resolution, 
which will require further large amounts of computing power in the future. Furthermore, 
we plan to produce so-called synthetic observations. These synthetic observations will 
allow us to directly compare our simulation results with real observations. Such compar¬ 
isons allow us to assess what can be inferred from observations - in particular how reliable 
parameters obtained from real observations are. 
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